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Abstract 

This paper proposes a centralized decision making framework at the macro base station (MBS) for device to 
device (D2D) communication underlaying a two-tier cellular network. We consider a D2D pair in the presence 
of an MBS and a femto access point, each serving a user, with quality of service constraints for all users. Our 
proposed solution encompasses mode selection (choosing between cellular or reuse or dedicated mode), resource 
allocation (in cellular and dedicated mode) and power control (in reuse mode) within a single framework. The 
framework prioritizes D2D dedicated mode if the D2D pair are close to each other and orthogonal resources are 
available. Otherwise, it allows D2D reuse mode if the D2D satisfies both the maximum distance and an additional 
interference criteria. For reuse mode, we present a geometric vertex search approach to solve the power allocation 
problem. We analytically prove the validity of this approach and show that it achieves near optimal performance. 
For cellular and dedicated modes, we show that frequency sharing maximizes sum rate and solve the resource 
allocation problem in closed form. Our simulations demonstrate the advantages of the proposed framework in 
terms of the performance gains achieved in D2D mode. 
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I. Introduction 


F uture fifth generation (5G) eellular networks are expeeted to be highly heterogeneous in arehi- 
teeture, with eoexistenee of maeroeells and femtoeells as well as deviee-to-deviee (D2D) eommu- 
nieations 0 - 0 . In partieular, femtoeells are of great importanee sinee they are predieted to generate up 
to 50% of the voiee ealls and up to 70% of the mobile data traffie in the near future Q. The two-tier 
eellular network arehiteeture, eomprising of a eentral maeroeell base station (MBS) and lieensed shorter 
range femtoeell aeeess points (FAPs), signifieantly improves the throughput for in indoor environments 
as well as the overall network speetrum and energy effieieneies 0, Q. However, eross-tier interferenee 
needs to be properly managed in sueh networks 0. 

Reeently, D2D eommunieations allowing direet eommunieation between nearby users has been envis¬ 
aged in 3GPP standards 0. The D2D users ean utilize unlieensed speetrum (out-of-band) or lieensed 
speetrum (in-band). Compared to out-of-band, in-band D2D ean provide more quality of serviee guar¬ 
antees 0 and is eonsidered in this paper. In in-band D2D, there are three modes of operation for the 
D2D users: (i) dedicated (or overlay) mode where D2D users are alloeated dedieated speetrum, (ii) reuse 
(or underlay) mode where D2D users reuse existing speetrum resourees and (iii) cellular mode where 
the D2D users are treated as normal eellular users and relay eommunieations through the MBS. From 
an operator perspeetive, determining the type of D2D operation during mode selection (assuming that 
neighbor diseovery has already been aehieved 0) is a erueial initial deeision by the network and an 
important researeh topie. In dedieated or eellular mode, the fundamental researeh ehallenge is resouree 
alloeation. In the reuse mode, the fundamental researeh ehallenge is interferenee management via effieient 
power eontrol. Overall, in order to provide operator managed quality of serviee guarantees, eentralized 
solutions whieh have low eomplexity are desirable. 

Literature review: Mode selection schemes have been proposed in the literature based on minimum 
distanee between the D2D transmitter (DTx) and D2D reeeiver (DRx) 0, biased D2D link quality and 
whether it is at least as good as the eellular uplink quality p0| or guard zones proteeting the MBS pT| 
or D2D users [12|. A limitation of the sehemes in 0, [10|, [12| is that they do not inherently proteet 
the D2D link from interferenee, while the seheme in pT| does not impose any restrietions on the D2D 
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distance; generally D2D eommunieation is envisaged as short range direet eommunieation. Also in [ 12[, 
the guard zone region surrounding D2D users is primarily used to determine whieh eellular users are 
allowed to reuse resourees alloeated to the D2D users, rather than speeifieally a mode seleetion eriterion. 
Mode switehing and mixed mode approaehes, where multiple modes are utilized at onee, are also studied 
in (Tg. 

Onee a mode is deeided, the network must address resource allocation to meet network requirements. 
In the reuse mode, power control is used to manage transmit powers and henee interferenee. Power eontrol 
is not guaranteed to provide closed form analytical solutions, but it has been shown that optimal solutions 


ean at least be found from searehing from a finite set [15|, although this elaim has only been made with 


two transmitting sourees in the system model. In |I6|, power optimization for one D2D transmitter and 
one cellular user transmitting during uplink was studied. Sinee there are two transmitters, the optimization 


is a simple two-dimensional problem. Power alloeation for maximizing sum rate was also studied in [ 17|, 
where the authors foeused on a binary power deeision, i.e., powers operate either at their maximums or 
minimums, and with no signal to interferenee plus noise ratio (SINK) guarantee for any user. The authors 
showed that binary power eontrol is optimal for two users, but is suboptimal for arbitrary number of users. 


Orthogonal resouree alloeation for D2D was studied in 1151 for both dedieated and eellular modes using 
the downlink (DL). In eaeh mode, time and frequeney alloeation was eonsidered, and for eaeh alloeation, 
greedy (uneonstrained) and rate eonstrained optimization was presented. However, in the rate constrained 
ease, only the eellular user has a minimum rate requirement, and thus the possibility exists for the cellular 
user to be alloeated all the resourees and leaving the D2D with none. Further, p3| only eonsidered a 
single-tier network in its system model. In a two-tier eellular network, a lieensed femtoeell ehanges the 
way resourees ean be alloeated, and in turn ehanges the maximization of the optimization objeetive in a 
non-trivial manner. Note that some papers use joint optimization [18|-[20| and/or game theory pT|-[25 1 
to solve resouree alloeation problems. In theory, joint optimization solutions eould be optimal, but their 
eomplexity often means approximations are required in praetiee. Further, in two-tier networks, different 
users may have different eonstraints or requirements whieh will further inerease the diffieulty of finding 
optimal solutions. Meanwhile, game theory has the advantage that it is a more distributive approaeh, but 
does not provide operator managed quality of serviee guarantees. 
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Both uplink (UL) and downlink (DL) spectrum resources can be used by in-band D2D. In the literature, 
there exists works which either use UL [ fT^ , [ fT^ , [ |2^ or DL ^n\ , [ [2^ , and also some which consider 


both p5| |, p9| . Generally, interference scenarios are less severe in the UL p2| , p0[ |, pT| . However, in 
this paper we assume DL resources are reused as this represents the worst case interference scenarios Q 

In summary, existing works on D2D communications have generally considered mode selection, resource 
allocation and power control sub-problems either separately or considered a subset of these problems 
for single and multi-tier cellular networks [(9|-p^, p3||-p9|, pT] |, p2| , [[2^-p9|. To the best of our 
knowledge, a centralized solution for mode selection, resource allocation and power control in D2D- 
enabled two-tier cellular networks is still an open problem. 

Paper Contributions: In this paper, we propose a base station assisted D2D decision making framework 
(cf. Fig. 2) that incorporates mode selection, resource allocation and power control in a two-tier cellular 
network. The MBS first decides if D2D dedicated mode is permissible or not based on the DTx-DRx 
separation distance being small enough and the availability of orthogonal resources. If not, an interference 
criteria is then applied to determine whether the D2D pair should enter reuse mode or remain in cellular 
mode. Resource and power allocation is then applied to maximize user sum rates. Compared to joint 
optimization methods, this multi-stage decision process can arrive at the correct mode and resource 
allocation in a much more straightforward fashion with less complexity. The major technical contributions 
of the paper are as follows: 

• We propose a mode selection method that prioritizes D2D dedicated mode if the D2D pair are close 
to each other and orthogonal resources are available, and otherwise allows reuse mode if the D2D 
pair satisfies a strict distance and interference criteria. We show that our proposed decision making 
framework allows more dedicated D2D users than conventional methods, and allows more correct 
decisions (i.e., higher rate) when resources are shared. 


• For the D2D reuse mode, we (non-trivially) extend the method described in [ 151 to three dimensions to 
solve the power allocation problem in a two-tier cellular network. In this process, we first analytically 
prove that (i) sum SINR is quasi-convex in any number of varying powers and (ii) sum rate has the 

'The methodologies developed in this paper can also be applied to reuse UL resources. We would only need to make a distinction between 
the two for cellular resource allocation since we assume half-duplex communications. 
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same derivative behaviour as sum SINK (and hence is almost quasi-convex) when one received 
power dominates in magnitude over others. Then using these results, we propose a simple approach 
of finding the corners or vertices of the power region to solve the power allocation problem, which 
achieves near-optimal performance as compared to exhaustive search. 

• For the cellular and D2D dedicated modes, we show that frequency allocation results in higher rates 
than arbitrary or time sharing resource block allocation. We solve the frequency allocation problem 
in a two-tier cellular network to maximize the sum rate, while meeting a minimum rate constraint for 
all the users to ensure fairness. We also present general resource allocation methods, where possible, 
for arbitrary number of users and transmitters. 

Paper Organization: This paper is organized as follows. Section presents the system model. The 
proposed framework and mode selection scheme is described in Section |in| The formulation and solution 


for power allocation problem in reuse mode is presented in Section IV The formulation and solution for 
resource allocation problem in dedicated and cellular modes is presented in Section |V| The results are 
presented in Section VI Finally conclusions are presented in Section |VII 


II. System Model 

We consider a single cell in a two-tier cellular network, as illustrated in Fig. Our system model is 
comprised of: (i) an MBS located at the center of the cell, which is serving a single cellular user equipment 
(CUE), (ii) an FAP serving a single femto user equipment (FUE), and (ii) a D2D pair comprising of a DTx 
and a DRx located close to each other. All the different user equipments (UEs), MBS and EAP are equipped 
with single omni-directional antennas. We assume that suitable inter-cell interference control mechanisms, 
such as fractional frequency reuse, are employed to avoid or manage inter-cell interference [j^. Hence, 
we study the single cell scenario. Although we study the simplified scenario, as illustrated in Pig. the 
proposed framework and resource allocation methods in this paper are applicable to the general scenario 
with multiple UEs and multiple PAPs, and will be discussed in their respective sections. A simple setup 
was also used in [|3^, with an included discussion on extending the model to more users. 


We assume that the MBS has perfect instantaneous channel state information (CSI) of all the links. 
This assumption has been widely used in the D2D literature pT]|, p^, [15|, |18| and allows benchmark 
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DTx goRx 


Fig. 1: System model comprising of a D2D pair, MBS, FAP, and its served users. Strong interferences to the DRx from the 
MBS and FAP are shown in red dashed lines. 

performance to be determined. The mode selection, resource allocation and power control is performed by 
the MBS in a centralized manner, based on the available perfect CSI. The transmit power of all transmitter 
nodes is denoted as Pt and the maximum transmit power is denoted as where t G {T,M,F} is 

the index for the transmitters, and T denotes DTx, M denotes MBS and F denotes FAP. The (minimum) 
rate at a receiver is denoted as (P™") TZr, while the corresponding (minimum) SINK under normalized 
resource allocation is denoted as ( 7 ™") 7 ,-, where r G {R, C, E} is the index for the receivers and R 
denotes DRx, C denotes CUE, and E denotes FUE. All the links are assumed to experience independent 
block fading. 

The instantaneous channel coefficients are composed of small scale fading and large scale path loss 
denoted as 

9t,r 

where n is the path loss exponent, ht^r is the small scale Rayleigh fading coefficients, which are assumed 
to be mutually independent and identically distributed (i.i.d.) complex Gaussian random variables with 
zero mean and unit variance and dt^r denotes the distance in meters between transmitter t G {M,T,F} 
and receiver r G {C, R, E}. For simplicity, we denote the distance between DTx and DRx as simply 
d. All links experience additive white Gaussian noise (AWGN) with power 

We use sum rate as our system performance metric with individual maximum power and minimum 
rate requirements. For clarity, due to the different nature of the D2D modes, we define each problem 
formulation in their respective sections. 
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Fig. 2; Proposed MBS assisted D2D decision making framework for mode selection, resource allocation and power control in 
D2D enabled two-tier cellular network. 


III. Proposed Framework and Mode selection 

We propose a base station assisted D2D decision making framework, as illustrated in Fig. to enable 
the MBS to deeide on the eorreet mode of D2D transmission, and determine the resouree parameters of 
whiehever mode is ehosen (power for reuse mode, frequeney resourees for dedieated and eellular modes) 
that will maximize sum rate subjeet to maximum transmit power and minimum reeeiver rate eonstraints. 
The main steps in this proeess are deseribed below: 

1) The MBS first deeides whether a potential D2D pair is elose enough for D2D eommunieations. 
Depending on the availability of orthogonal resourees and potential interferenee, dedieated or reuse 
mode is ehosen. Otherwise, the pair remain in eellular mode. 

2) If the reuse mode is ehosen, then the MBS instruets the CUE, DTx and FAP to eontrol their transmit 
powers to guarantee quality of serviee to all reeeivers. This is done aeeording to the approaeh 


proposed in Seetion IV 


3) If the eellular mode or the dedieated mode is ehosen, then the MBS alloeates resourees to the CUE, 
EAP and D2D UEs. Sinee interferenee is not present, all transmitters ean use the maximum transmit 
power. In the dedieated mode, we assume D2D UEs use the DE resourees. In the eellular mode. 
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we assume that both UL and DL resources are used by D2D UEs since the D2D communication is 
being relayed by the MBS. 


A. Mode Selection 

In order to allow as many potential D2D dedicated users as possible, our decision making framework 
firstly allows a potential D2D pair to enter dedicated mode if they are close enough and orthogonal 
resources are available. Interference is not considered in this case as using orthogonal resources eliminates 
interference. If orthogonal resources are not available, the decision to enter reuse mode is then made based 
on the potential interference. A potential D2D pair can only enter reuse mode if both of the following 
criteria are satisfied: 

1) The DRx must be located outside an interference region such that the potential interference is lower 
than a threshold. Since actual powers are yet to be determined, we assume maximum transmit 
powers. 

2) The distance d between the DTx and DRx must be less than a threshold (should be satisfied from 
initial step). 

To determine the distance and interference thresholds, we recognize that in the cellular mode, the CUE 
and EUE should always experience better rates compared to D2D due to less interference. Therefore, 
mode selection is equivalent to finding under what conditions the DRx SINR in cellular mode is better 
than in D2D mode for the D2D users, i.e., 

T 9t,m 9m, r \ y -fr S't.r 


mm 


( 2 ) 


P 'max„ I _2 ’ pmax„ i _2 / — pmax„ i pmax„ i _2 

F 9f,m + ct A o j ^jm.r + -fp 9f,r + cr 

where the min(-, ■) denotes the minimum operator and is used since the rate of the cellular two-hop link 
is limited by the minimum of the UE and the DE. 

Suppose we consider a scenario where a D2D pair is close to each other, but located within a high 
interference region. Erom Q, if the interference is greater than a certain threshold 


Pm''9m,r + Pr9F,R > Pr9T,R X min 


pmax^ I —2 pmax^ i — 

-r F + CT Jr'p 


- CF , 


( 3 ) 


pmax „ ’ pmax „ 

TT 9T,M 9m,r 

using D2D mode will be an incorrect decision as it will lead to a lower rate. Note that this threshold 
value is a conservative estimate as it does not consider the overall system improvement from the CUE or 
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FUE. When the eorreet mode seleetion deeision is made, the rate aehieved by the D2D pair, and also the 
overall system, will be greater. 

A similar argument ean be made for a seeond seenario where the DTx and DRx are loeated outside a 
high-interferenee region, but are far apart. Rearranging Q to solve for the D2D separation distanee, we 
find that D2D mode would be an ineorreet deeision if the DRx is outside the interferenee region, but the 
D2D pair is separated by a distanee of 


^adaptive ^ 


'-T.R 


pmax 

-^D 


{P^^^9m,r + Pr9F,R + ^ 2 ) 


X 


' mm 


( pr"9F,M + 

V pr9T,M 


Pr"9F,R + a^ \ 
Pm''9m,r )' 


( 4 ) 


where n is the path loss index. 

However, when the DRx is elose to an interferenee souree, Q may provide an unneeessarily small 
threshold, and therefore limit the number of D2D pairs. Thus, in our framework we ehose the maximum 
threshold between Q and a predetermined value ^constant, i-e.. 


d ^ maX-[(iconstant! *^adaptive} 


( 5 ) 


The benefits of this approaeh will be illustrated using simulation results in Seetion VI-A 


IV. Power Allocation in Reuse Mode 

In this seetion, we solve the overall sum throughput optimization problem in the reuse mode. In reuse 
mode, the problem reduees to finding the optimal powers that ean maximize the sum throughput objeetive 


while meeting individual minimum rate requirements. We extend the method in [15| for the ease of two 
transmitters to three transmitters to solve the power alloeation problem in a two-tier eellular network. We 
present a geometrie representation of the problem for the ease of three transmitters (i.e., DTx, MBS and 
FAP) and present a near-optima^ solution approaehing that of exhaustive seareh. 

^ We use the expression “near optimal” to describe the closeness of the solution to the optimal solution, rather than to define it as a 


specific solution or class of solutions. 
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A. Problem Formulation 


Our overall system aim is to maximize the sum rate with individual transmit power and reeeiver rate 
eonstraints. We ean formulate the optimization problem as follows: 

-Pt5't,r ^ I 1 A I -Pm5'm,c \ 

( 6 ) 


max < 7^ = log 2 1 + — - ^ + log 2 1 + , 

Pt,Pm,Pf [ \ PMgM,R + PFgF,R + J \ -Pt5'T,C + -< f5'F,C + C’’ 


+ log 2 ( 1 + 


-PfS'f,! 


-PtS't.e + -Pm(?m,e + 


sueh that 


Pt < Pr"", t e {M, T, F} 


-Pt5't,r 

PMgM,R + -PfI/F.R + 
-Pm5'm,c 

-Prfi'T.c + -Pf5'f,c + 
-PfS'f.e 

-Prfi'T.E + PMgM,E + cr^ 


>7r 


>7^ 


>7r 


(7a) 

(7b) 

(7c) 

(7d) 


where ( |7a| ) represents the maximum power constraints for each transmitter, while fTb])—(|7dl) are minimum 


SINK requirements. Note that for reuse mode, since all resources are shared and allocation is not 


considered, a minimum rate constraint is equivalent to a minimum SINK constraint. 


B. Geometric Representation 

We adopt a geometric approach to determine the optimal powers. To graphically represent the admissible 
powers, we first set orthogonal axes to be the powers. Next, setting constraints (|7b|)—(|7d|) to equality and 
rearranging, we obtain 

/t — 5'T,R-Pr — 7 r™5'm,r-Pm — 7r'"5'f,r-Pf — 7™"*^^ = 0) (8a) 

/m — — 7c'"5't,c-Pt + 5'm,c-Pm — 7c'"5'e,c-Pf — 7c'"cr^ = 0) (8b) 

/f — —7e'"5't,e7^t — 7e™5'm,e7^m + fi'F,E7^F — 7™"*^^ = 0, (8c) 

which represent planes in 3-dimensional space. The planes themselves represent the relationship between 
each node’s power and the SINK thresholds. Each plane focuses on one threshold, and thus we refer to 


(|8 a|)— as the D2D, MBS, and FAP planes respectively. Each plane intersects with its respective axis at 
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their respective minimum powers P™". Note that while the thresholds are stated in terms of the receiving 
node of that link, the powers are of the transmitting node. 


We can plot ([8^—([8^ using their inequalities to obtain a 3-dimensional upper right comer region 
within a cube|^ the faces of which represent the maximum individual power constraints. The top right 
comer of this cube has the maximum power coordinates (P™, P™“). 

The smallest possible transmit powers, P™“, that satisfy each users’ SINK requirement can be calculated 
from fTb])—fTd]) when there is no interference from the other transmissions. Therefore, the range of 
admissible powers is 


min 2 

^min _ /R ^ ^ pmax 

9t,r 

min 2 

min / C ^ ^ TD ^ omax 

S EM S Pvi 


pmm _ 

-'^M — 




min 2 

^min _ /E ^ ^ P pmax 

■'F — F — -'F' 

fi'F,E 


(9a) 

(9b) 

(9c) 


Meanwhile, the minimum powers that jointly satisfy the individual user rate constraints can be found by 


simultaneously solving (|8a|)—(|8cl) using standard methods such as Cramer’s rule. Note that these powers 
will not maximize sum rate. 


We assume that the coefficient matrix formed from ([8a|)—([8^ is full rank, i.e., the three planes intersect 
at a point Q, whose coordinates are all positive values since they represent physical transmission powers. 
We find that the reuse mode is a viable option only if each signal strength is relatively large compared 
to the interference, making it easier to satisfy SINK constraints. This conclusion is consistent with others 


in the literature p4| |. 

The admissible power region is formed by the intersection of the three planes in 3-dimensional space, 
and is bounded by these three planes and the three faces of the cube. The optimal powers lie within this 
power region. In order to avoid an exhaustive search, which would be computationally expensive, we 
propose a near optimal solution which reduces the process to testing and selecting the optimal powers 
from a finite set. 


^Strictly speaking, the region is a rectangular prism, but for conciseness we will use ‘cube’ to describe this region. 
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(a) All thresholds are satished. 



(b) Two thresholds are satished. 



(c) One threshold is satished. 


Fig. 3: Illustration of different types of power regions. 


C. Proposed Solution - Vertex Search 

In this paper, we adopt the simple approach of finding the comers or vertices of the power region to 
test for the optimal powers. This approach relies on the following two mathematical conditions: 

1) The optimal powers cannot lie in the interior of the power region, and must be on a boundary. 

2) The objective function is quasi-convex on a boundary, ensuring that the maximum values are at the 
endpoints/vertices. 

The first condition was in fact proved in pT| , and thus it is known that at least one of the powers is 
at its maximum when maximizing sum rate. However, this only states that the optimal solutions exist on 
the boundary of the power region, which includes vertices as well as higher dimensional edges and faces 


that contain an infinite number of points. Thus, this conclusion from [17| alone is not sufficient to obtain 


the finite set of points which will give the optimal solution. For two transmitters, it has been proven that 


the optimal power lies on the corners or vertices of the power region [15|, [16|, a fact that relies on the 
convexity of the sum rate function for two powers. However, it is well known that in general, the sum rate 
expression in @ is non-convex with respect to arbitrary combinations of varying powers. Consequently, 
for arbitrary number of transmitters, the optimal powers may not necessarily lie on the vertices of the 
power region, leading to a possibly infinite set of points to test. 
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To prove the second condition and justify searching the vertices to maximize sum rate for arbitrary 
number of powers, we present the following two propositions. 

Proposition 1. Sum SINK is a quasi-convex function for any combination of varying powers. Hence, it is 
also jointly quasi-convex in all powers. 

Proof: See Appendix ■ 

Remark 1. Since sum SINK is a quasi-convex function, the powers maximizing it will lie on the one of 
the vertices of the power region. 

Proposition 2. When one receive power dominates, global maxima and minima for sum rate and sum 
SINK will occur at the same locations. 

Proof: We prove in Appendix that when one receive power dominates, e.g., an order of magnitude 
larger than others, sum SINK in ( [2^ and the inner log term in pT] ) have the same asymptotic derivatives, 
meaning that the two functions will ‘follow’ each other more and more closely the larger the dominant 
power is. Since logarithm is a monotonic function and does not change the locations of local maxima or 
minima, this implies that the same powers that maximize sum SINK will also maximize sum rate. ■ 

Remark 2. Since global maxima of sum SINK will be at the vertices of the power region, the same vertices 
will also give near-optimal solutions for sum rate. 

Note that approximations such as reformulating the objective function as a geometric program (GP) 
p5| can be used to solve However, we show in the results section that our proposed simple approach 
yields near optimal solutions quite close to those obtained using exhaustive search and GP, but does not 
require an iterative approach. 

In the next subsection, we present a systematic way of obtaining the coordinates of the vertices of the 
power region by solving relevant sets of SINK equations. 

D. Vertices of the Power Region 

We identify all the vertices of the power region for the different interference scenarios, from which one 
of the points will give us the suboptimal powers that maximize the sum rate. All the vertex points are 
summarized in Table [ij The notation {Pa, means solve for powers Pa and Pf, using simultaneous 
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TABLE I: Finite set of vertices (suboptimal powers) for reuse mode. 


Type 

Condition 

Number 

of points 

Set of vertices (suboptimal powers) 











Face point 


All 

9 

({Pt, Pf} y, PS“), ({Pt, Pf} y, PS“), 





(pr\ {J’m, PF}y^^’^“>),(pry {^’m, PF}y^^'^^>), 





(P^?“,{J’M,PF}y^“’^^>) 

Edge point - 

All thresholds 

7r ^ 7r 5 7c ^ 7c 

4 

(P^|{/t}^ pmax^ pmax)^ (pmax^ Pm 1 ^^“ > , Pr^, 

satisfied 


and 7^5 > 7 g--' 


(p-f^L ps“, pfI^^'^^), (Pry -Pry ^ry 



7e < 7b ” 

4 

{Pt\^m,p^^\ pry, (PT^^-^y pry ^^ry 





(Pry PM^^^y pry, (^ry PM^^'^y pry 

Edge point - 

Two thresholds 


7r < 7r " 

4 

(Pry PM^^^y pry, [p^\ PMy^^y pry 

satisfied 




(pry pry PF^-^'^y, (pry ps^y PF^^'^y 



7c < 7c‘" 

4 

(PTy^T>, PiS“, Pf“), (Pxy^^y pry j^ry 





(Pry pry Pry^^y, (pry ^^ry Ppy^^y 

Edge point - 

7i5 

< 7b " and ic < 7c " 

2 

(PTy^“>, pry ^ry, (Pxy-^^y pry p'ry 





One threshold 

satisfied 

7r 

< 7r and 7 ^ < 7 “" 

2 

(p)?^, Pm y y pry, (p'x“, em y y pry 

7r 

< 7 S‘" and 7 t < 7 ^“ 

2 

(Pry pry Ppy^^y, (pry ^^ry Pxy^^y 





IEEE TRANSACTIONS ON COMMUNICATIONS SUBMISSION 


14 


equations fa and fb with the other power maximized, where a,b G {T,M,F}. How these are obtained is 
explained in detail below: 

Face points with one power maximized: There exists vertiees that lie on a faee of the eube and are 


formed from the interseetion of two planes, e.g., point F in Fig. 3a There are nine sueh vertiees (three 
faces with three ways of choosing two intersecting planes for each face). These vertices can be found by 
solving two plane equations simultaneously with the power corresponding to the third face maximized. In 
general, it is difficult to identify exactly which of these nine points may be optimal for a given interference 
scenario. Thus, we need to test all nine vertices. 

Edge points with maximum powers satisfying all thresholds: Consider the case where the three planes 
are orthogonal, as shown in Fig. In this case, the power region includes the top corner of the cube, 
where all three powers are maximized, and three other corner points where the planes intersect the edges of 
the cube, which we shall label as edge points. Since the top comer lies in the power region, this indicates 
that when all powers are maximized, all three SINRs 7^ are greater than their minimum thresholds 7™". 
For the rest of this section, we denote 7 ^ as the SINK for each node when all powers are at their maximum. 
There are four such points, as summarized in Table I. Note that the same SINK scenario can occur even 
when the planes are not perpendicularj^ The distinctive feature of this scenario is that the top comer is 
within the region spanned by the planes, and that each plane only intersects one of the maximum power 
edges of the cube. 

Edge points with maximum powers satisfying two thresholds: To visualize this scenario, imagine 
tilting the planes pivoted at Q to form new power regions. For instance, if we tilt only the FAP plane 
upwards, it will eventually pass through the top corner and intersect the other two top edges. These two 


additional edge points {B and D in Fig. add to the existing two edge points {A and C) to give a total 
of four edge points on the cube’s edges. Since the top comer point will now be below the FAP plane, 
this means that 7 ^ < 7 ^'". Similar arguments can be made for the other two planes, giving us three cases 
where there are a total of four corner points in the power region, each case corresponding to one 7' that 
is less than its respective threshold. 

Edge points with maximum powers satisfying one threshold: For scenarios where two 7' fail to reach 

''in fact, perpendicular planes which each only intersect one axis corresponds to an interference-free scenario. 
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their thresholds and only one is met, the two planes will be tilted sueh that the eorner point lies outside 
both their feasible regions, as shown in Fig. where the FAP and MBS planes lie above and to the left 
of the top eorner respeetively. In these eases, the feasible region will interseet one of the three eorner 
edges at two points. Fig. [^illustrates the maximum MBS power edge being interseeted at two points A 
and B by the D2D and FAP planes respeetively. Note that the MBS plane also interseets the same edge, 
but that point of interseetion is outside the power region. Thus, we get two points, eaeh eorresponding to 
a set of eonditions. 


V. Resource allocation in dedicated and cellular modes 


If mode seleetion deeides that the D2D pair ean transmit using either dedieated or eellular mode, time 
and/or frequeney resourees must be alloeated. We make the following assumptions for resouree sharing 
in both dedieated D2D and eellular mode: (i) sinee eellular frequeneies are used, there is a minimum rate 
guarantee for eaeh user, ineluding the DRx, (ii) there are enough resourees to meet all users’ minimum 
rate requirements, and (iii) at any one time, one transmitter ean only operate in either uplink or downlink. 


i.e., half duplex. Consideration of full duplex transmitters [36| is outside the seope of this work. 


A. Problem Formulation 

Sinee there is no interferenee in both dedieated and eellular modes and all powers ean be maximized, 
the SINR at eaeh reeeiver is the same as the signal-to-noise ratio (SNR) at that reeeiver, given a^ 


7 c = 


9m,cPi 


max ^ pi 

M _ 9T,RFr 

2 -’ - 


max ^ pi 

T _ 9f,eFi 

2 —’ - 


max 

F 


( 10 ) 


We formulate a general optimization for a Long Term Evolution (LTE)-like resouree grid with distinet 
resouree bloeks as follows 


A* 


maximize 


subjeet to 




Bi6f 


A* 




( 11 ) 

( 12 ) 


^With slight abuse of notation but for the sake of simplicity, we use the symbol 7 r for SNR, where as in Section II we denoted minimum 
SINR at a receiver as 7 ^^^, where r € {R, C, E} is the index for the receivers. 
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where 5* is the number of resource blocks for user r at the ith. time interval, 5^ and 5l are the (constant) 
fractions representing the portion of each frequency and time block compared to the total grid respectively, 
and iV* is the total number of time intervals for the resource grid. The total number of blocks allocated 
to user r is therefore Note that we divide the SNR by the frequency portion as we define 

with respect to the entire bandwidth, resulting in equal noise power density [13]. 

The general formulation is difficult to solve, and in practice requires numerical methods. In order to 
gain insight into generalizations for arbitrary number of users and to obtain closed form solutions, we can 
show that for a given B^. = M number of resource blocks and assuming that each block is no more 
preferable to any other, allocating resources across frequency will produce higher rates than allocating 
across time or in a random manner. For ease of proof and without loss of generality, we assume high 
SNR such that log 2 (l + SNR) ^ log 2 (S'iVi?). Using log 2 (A) + log 2 (i?) = \og 2 {AB) and the arithmetic- 
geometric mean inequality, which states that the maximum of a product of terms with a sum constraint 


occurs when all terms are equal, we find that the maximum of 



Aft 

subject to E Bl=Rf 

i 


(13) 

(14) 


occurs when all 5* are equal. In other words, using equal bandwidth allocation across all time intervals 
for each user will provide the largest rates. Thus, although we can generalize resource allocation to be 
compatible with arbitrary resource block allocations, frequency allocation will give higher rates compared 
to other approaches for a given number of resources blocks. 

Since frequency allocation can be solved in closed form, we analyze frequency allocation formulated 
as follows: 


maximize 

Xr 

subject to 


5^a;,log2 + 

a;,log2 + 


(15a) 

(15b) 


where the factor Xr is a function of the resource portions 0 < a, a', 13,13' < 1. For ease of analysis, we 
study only the case where an exact bandwidth is allocated across all time intervals. 
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(a) Frequency sharing in dedicated mode. 


(b) Frequency sharing in cellular mode. 


Fig. 4: Illustration of resource allocation in dedicated and cellular modes. 


B. Frequency sharing in dedicated D2D mode 


For the allocation structure illustrated in Fig. 4a we want to maximize 


112 = a log2 (^1 + + a' log2 (^1 + +{l-a-a') log2 (^l + 


7e 


(16) 


1 — a — a' 

1) Unconstrained: With no minimum rate constraints, we differentiate ( [T^ with respect to a and a', 
and simultaneously solve for = 0 and = 0, which gives us the solutions 


a = 


7c 


7c + 7 r + 7e ’ 


a = 


7r 


7c + 7 r + 7e 

Substituting the above into ([T^ and simplifying, the optimal sum rate is 


(17a) 

(17b) 


'nT = log 2 (l + 7c + 7 r + 7e)- 


(18) 


2) Constrained: To meet each user’s minimum rate requirement, we require the solution to 

«log2(l + ^) =7^^ (19) 

for each user. The solution can be written in terms of the Lambert W-function (see Appendix |^, but 
there is no analytical solution that can be expressed using elementary functions. A simple numerical line 
search along 0 < a < 1 can be used to find an optimal solution. 
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C. Frequency sharing in cellular D2D mode 


In frequency sharing cellular mode, because there is only one MBS transmitter, the D2D UL and CUE 
UL must occur at the same time, with D2D DL occurring immediately afterwards. The FUE can be 
allocated subbands at any time as it is served by a separate transmitter. Therefore, the allocation scheme 


follows the partitions as shown in Fig. 4b 


In frequency sharing in cellular mode, the sum rate to be optimized is 


7m,UL 
a 


+ mm 

7e 


7^4 = Ol(3 log2 (l + 

+ {l-a-a) logs + 7-. - — - 

\ [1 — a — a') J 


'/3 logs (l + ’ (“ + " logs (^1 + 


7r,dl 
a + a' 


( 20 ) 


(1 — a — a') / ’ 

m pmax 

where /7+/7' = and 7r,uL = — is the SNR at the MBS during CUE UE with the CUE transmitting 

at its maximum power 

1) Unconstrained: We define uplink and downlink rates as 

TR.DLA ( 21 a) 

(21b) 


T^ul — logs (^1 + 
T^dl = logs ( 1 + 


a' 

7r,dl 
a + a' ^ 

To simplify ( |20] ) into an expression involving only a and a', we note that the maximum sum rate occurs 
when a'/37l\ji, = (a + with the solution given by 

T^dl 


/5 = 


OL+a 


yT^UL + 7^: 


( 22 ) 


DL 


Substituting the above, the rate expression for D2D is 

ci'T^ulT^dl 


'Tldia.) = 


jT^ul + 7?.dl 


(23) 


a-\-a‘ 


Therefore, we can simplify ( [20| ) to 

«7^DLl0g2(l+2^^: 


— 


a+Q' 


jT^ul + 7^: 


+ TZd{a) + {I - a - a') logs ( 1 + 


DL 


7e 


(1 — a — a') 


(24) 


A numerical search for 0 < a + a' < 1 can be performed to maximize (24). 


2) Constrained: In this scenario, we desire to maximize ( |24| ) under a minimum rate constraint for each 
user. Again, a numerical search can be performed to find the maximum. 

^Setting /?' = 0 would be equivalent to having two cellular users, and the solution would be the same as that in Section V-B. 
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Fig. 5: Percentage of potential D2D pairs entering dedicated mode. Predetermined threshold is better when interference is 
large, while adaptive threshold is better when interference is small. 


VI. Results and Discussion 


In this section, we present simulation results to illustrate the benefits of using our decision making 
framework over conventional cellular transmission for a potential D2D pair. Unless stated otherwise, 
simulation parameters presented in Table are used, which are similar to those adopted in |18|. We use 
(x, y) coordinates in meters to describe node locations. 


A. Mode Selection 

We first show the advantages of using Q compared to using a constant and adaptive distance threshold 
only. Setting ^constant = 50m and assuming orthogonal resources are available 50% of the time. Fig. shows 
that picking the largest threshold between the predetermined and calculated gives the highest percentages 
of users entering dedicated mode. When interference to the DRx is large, choosing a predetermined 
distance threshold is more beneficial. When the DRx is farther from an interference source, an adaptive 
threshold is the better choice as larger D2D separation distances can be tolerated. It is evident that the 
proposed method captures the best features of the other two, and in fact slightly outperforms the best of 
both at every location tested. 

When orthogonal resources are not available. Fig. [^plots the D2D rate gain versus the distance between 
the DTx and DRx, d. The D2D rate gain refers to the ratio between the D2D rate and the cellular rate. 
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TABLE II: Values of Simulation Parameters 


Parameter 

Value 

Bandwidth 

20 MHz 

Noise speetral density 

-174 dBm/Hz 

Max MBS transmit 

= 43 dBm 

power 


Max FAP transmit 

= 21 dBm 

power 


Max DTx transmit 

= 23 dBm 

power 


DTx eoordinates 

Varying along x = y 

DRx eoordinates 

Varying along x = y 

MBS eoordinates 

(0, 0) 

CUE eoordinates 

(500, 0) 

FAP eoordinates 

(100, 200) 

FUE eoordinates 

(110, 200) 

DTx to DRx pathless 

28 + 401ogio(d) (dB) 

MBS to CUE pathless 

15.3 + 37.61ogio(dM,c) 


(dB) 

FAP to FUE pathloss 

38.5 + 201ogio(dF,£;) 


(dB) 

CUE minimum SINR 

7 ™" = 0 dB 

FUE minimum SINR 

= 7 dB 

DRx minimum SINR 

7 ™" = 3 dB 


both under the same interferenee eonditions. We ean see that the D2D gain deereases when the D2D pair 
beeome farther apart and also when the DRx is eloser to the MBS. This is in line with the diseussion in 
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(a) D2D rate gain versus the distance between the DTx and 
DRx, d for different MBS-DRx distance, (Im.r- The shaded 
area below D2D rate gain of 1 represents the region where 
selecting D2D mode would be an incorrect decision. 



(b) D2D rate versus the distance between the MBS and DRx, 
for mode selection using distance only criterion and 
two stage criteria. 


Fig. 6; Performance of proposed D2D mode selection criteria. 


Section III since: (i) when the DRx is closer to the MBS (the largest interference source), using cellular 
mode should provide higher rates than an incorrect D2D mode decision since there would have been 
more interference, and (ii) when the D2D pair separation distance increases, cellular mode should provide 
higher rates since D2D mode would be weaker with increasing separation distance under constant transmit 


power. In Fig. 6a the D2D separation distance at which each curve intersects the boundary of this region 
can be calculated using Q. Our calculated and simulated values were found to be in close agreement. 
For example, the calculated separation distance for dM,R = 600 m is 75.9 m, while the simulations give 
a value of 71 m. 


Fig. 6b shows the actual rates experienced by a DRx when using the proposed mode selection method 
satisfying @ and when using just the D2D minimum distance criterion with d = 50 m. If D2D mode 
is always allowed for d = 50 m, the DRx can experience a smaller rate due to its close proximity to an 
MBS (or other large interference source), while using our proposed method will avoid such instances. 


It is important to note that our results in this subsection do not suggest that cellular mode is superior to 
D2D mode. Rather, our results highlight that under some conditions, using a single criterion to determine 
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mode selection can lead to an incorrect decision. We will show in Section VI-C that if D2D is operating 
in dedicated mode, it does in fact outperform regular cellular mode. 


B. Reuse Mode 

Fig. 1^ plots the sum rate in reuse mode versus varying MBS-DRx distance, dn.R, comparing the near- 
optimal powers found using the proposed approach with those obtained from geometric programming 
(GP) [ |M| and exhaustive search. We believe that using GP and exhaustive search serve as sufficient 
benchmarks - GP is one of the most common numerical approaches to finding near optimal solutions 
for the power control problem, while exhaustive search with a sufficiently fine step size confirms the 
optimality. 

Our results show that for the considered parameters, our proposed method of searching the vertices of 
the power region and using the one that gives the maximum sum rate is comparable to optimal solutions. 
However, our method requires far fewer calculations than exhaustive search and GP, the latter of which 
relies on successive approximations with no prior indication on how many iterations are required. For 
example, using an Intel iV 3.2 GHz CPU with 16 GB RAM, for dn.R = 1000 m in Fig. 7 GP took 
up to 44.5 seconds to calculate a solution, exhaustive search took 49.8 seconds, while our vertex search 
took only 1.2 seconds, i.e., an improvement of around 40 times over both benchmarks. Thus, GP can be 
unreliable in determining a suboptimal solution in sufficient time, while our vertex search approach will 
always return a suboptimal solution if the problem is feasible for small numbers of reuse powers. A further 
advantage of vertex search is that it always takes approximately the same time to calculate a solution for 
each realization, while the run time and accuracy of GP heavily depends on stoppage parameters. 


C. Dedicated and Cellular Modes 

Fig. shows the sum rate gain, i.e., sum rate in dedicated mode divided by sum rate in cellular mode, 
under minimum rate requirements for each user. It is clear that dedicated D2D mode provides a greater 
sum rate when the D2D separation distance is small, and/or when the distance between the MBS and 
DRx is large. 

It must be noted that the unconstrained dedicated and cellular modes offered similar sum rates under 
the simulation parameters, and thus their results are not shown. It is clear however that unconstrained 
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Fig. 7: Sum rate in reuse mode with transmit powers deter- Fig. 8: Sum rate in reuse mode with transmit powers deter¬ 
mined using proposed near-optimal vertex search approach, mined using proposed near-optimal vertex search approach, 
geometric programming and exhaustive search. geometric programming and exhaustive search. 


dedicated sum rates will never be lower than their cellular counterparts since the D2D option is intended 
to improve overall system performance, and will not degrade the best performing user. Thus, we can 
conclude that D2D mode is more advantageous when users have individual rate constraints. 

D. Scability Discussion 

Although we have presented our methodology using a simple system model, we can analyze the 
scalability with respect to increasing number of base stations and users. For our mode selection framework, 
the number of decisions scales linearly with the number of potential D2D pairs, and not the total number 
of users or base stations. 

For reuse mode, we have presented three transmit powers to be optimized, leading to a 3-dimensional 
problem. Increasing the number of transmit powers increases the dimensionality of the problem, while 
increasing the number of users increases the number of planes, and further restricts the size of the power 
region. For N powers, the power constraints form an iV-dimensional hypercube, while the minimum SINK 
constraints further bound the region to form an A^-dimensional polytope. Depending on which scenario 
the network is in (i.e., number of thresholds satisfied by max powers), the complexity of vertex search 
could increase exponentially at worst (e.g., in Fig. 3a), and linearly at best (Fig. 3c). However, if a small 
number of users share the same resource, since our vertex search avoids iterations, it can still be a more 





IEEE TRANSACTIONS ON COMMUNICATIONS SUBMISSION 


24 


effective solution. Exact expressions for the vertices for iV-dimensions is an interesting topic for future 
research. 

For frequency sharing in dedicated mode (illustrated in Fig. 4a), we prove in Appendix that the 
unconstrained case has a general solution for any number of transmitters and partitions. The general 


constrained case also has a solution given by solving ( [191 ) for o^oh user. With increasing number of 
nodes (base station or user), the complexity would scale linearly as each additional node would require 
one additional equation to solve for (from differentiating ( [T^ or solving ([T9j)). Frequency sharing in 
cellular mode (illustrated in Fig. |4^ can have various allocation structures due to the simultaneous uplink 
condition, and thus cannot be generalized in the present manner. 


VII. Conclusion 

We have presented a comprehensive mode selection, power control and resource allocation framework 
for D2D communication underlaying a two-tier cellular network. Our proposed mode selection scheme 
allows D2D communications under stricter conditions, leading to more correct decision making and a 
higher rate of allowing dedicated mode. We have also proposed a geometric approach to determine 
near-optimal powers for power allocation in reuse mode with faster computational time than benchmark 
methods, and provided closed-form resource allocations for orthogonal D2D mode for any number of 
users. 

There are numerous interesting additional features and directions for future research in our work. Energy 
efficiency could be used instead of sum rate as an objective, which would be particularly relevant for 
uplink scenarios, while imperfect CSI is a practical issue that can also be considered. 


Appendix A 

Proof of Proposition d] 


From [371, a differentiable function is quasiconvex if and only if 


fivi, ■ • • ,2/n) < fixi, V/(a;i,. . .,Xn)'^{yi -Xi,...,yn- Xn)'^ < 0 


(25) 


Although we can apply this to the sum SINK function directly, we note that since the addition of bounds 
and differentiation are preserved under addition, it is sufficient to show that each SINK is quasiconvex in 
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order to show that sum SINK is quasiconve?([^ Further, we ean ignore the noise eonstant in the denominator 
as it does not ehange the eonvexity behaviour or shape of eaeh fraetion. 

Consider the generie definition of sum SINK 


N 


5 = E 


P, 


(26) 


2 = 1 


where Oj = + cr'^- For N varying powers, if the numerator of an SINK fraetion is eonstant, i.e., 

a power that is not varying, then 

k 

(27) 


Pi + . . . + PlSf 

for varying powers Pi,, P/v and eonstant k is elearly quasieonvex as it follows a hyperbolie shape. If 
the numerator is a varying power, i.e., 

(28) 


P2 + ... + Pat-i 


then using the seeond inequality in (25) we find that for two sets of powers {Pi,..., Pat} and {P(,..., P(^}, 


VP(Pl, . . . , P.)^(P{ - Pi, ..., P,; - Pnf = 


1 


P2 + ...+Pn 


-Pi 

(P2 + ...+Piv)2 


/ P( - Pi ^ 


\Pn — Pn J 


P[ - Pi 


N 




P' - P 


P2 + . . . + P/v (P2 + . . . + P/v)^ 

N 


< 0 , 


(Pi - p^m +... + fn) - -p. J2^p; - Pi) < 0 , 


i=2 


P{(P2 + ...+P/v) <Pl(P^ + ...+P^), 
Pj ^ Pi 

P2 + • • • + P'n P2 + ... + Pn 


(29) 


This is the first inequality in (25) when P/ = i/i. Pi = xi. Thus, for any eombination of varying powers. 


we find that sum SINK is quasieonvex. 

^In general, the sum of quasi-convex functions may not be quasi-convex. 
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Appendix B 

Proof of Proposition [2] 


Suppose we differentiate S in ( [2^ with respeet to the most dominant power Pi'. 

dS _ 1 _^Pj 

dPi tti ^ a? ’ 


(30) 


If Pi was a dominant power, we observe that the derivative will approach l/a* since Oj —)■ oo as P* —)■ oo. 
Similarly, for the generic definition of sum rate 

P \ . -rV / . P 


N 


P = X] log2 1 + 


2=1 


N / 

iog2 n ^+ 

i=i \ 


Pj + ^V *=i V 

if we expand out the brackets in pT) ignoring the logarithm to obtain 

^ / p\ ^ p. 

of SINRs two at a time) 

7=1 k ^ 


(31) 


N 


(Products of SINRs three at a time) + ... + 


P 


(32) 


2=1 


and differentiate with respect to p, we find that all the derivatives of the products of SINRs will contain 


aj in the denominator, and will approach 0 as p —)■ cxd. Thus, both sum SINR and (32) have the same 


asymptotic gradient of l/a* when one power dominates. Note that if we differentiate with respect to Pi, 
but p was not the dominate power, both expressions will instead approach —Pj/a'j if Pj was the dominant 
power. 


Appendix C 

Closed form solution for constrained frequency sharing in dedicated D2D mode 


To solve ( [T9l ), we need to manipulate ( fTO] ) to a form where we can use the Lambert W function. Firstly, 
we can rearrange and then exponentiate ([T^ to get 


a 


(33) 


Next, we need to introduce additional terms such that the exponent contains the left hand side, i.e.. 




a 


( 34 ) 
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Moving the exponential to the left hand side gives 


7^““ln2 


In (1 + ^ ) e 
a 


_ 77 ™'“ In 2 


(1+^) = 


In 2 


2-7^^“/7r_ 


'yr ^ OL J 

We can now apply the Lambert W function since the exponential is in the form Ae^: 

7^^“ln2, 


7^““ln2, . 7^ 


7r 


In (1 + ) = w 

a 


‘_2-nr'^hr 


It 


Rearranging for a gives the solution 


a = 


- 7 ^ 7 ^““ In 2 


In 2 + 7 ^W ( 

r " \ "fr 


To ensure a real solution, we use the —1 branch of the Lambert W function. 


(35) 


(36) 


(37) 


Appendix D 

General solution for unconstrained frequency sharing in dedicated mode 


For N transmitters, and hence N partitions, sum rate is 


N 




(38) 


where ^ = 1 is the partition fraction and 7 * is the SNR of each receiver. 

In order to greedily maximize R, we need to simultaneously solve the partial derivatives with respect 


to each at, i.e. |^ = 0. This will give the relations 


doi 


li = JA 

oti (y.j^ 


(39) 


for 7 , fc = 1,..., A^. Setting k = m, and noting that = 1 — J2k=i we can rearrange (|39[) to obtain 


CXi 




C77fc \ li li 


Oti 


k=l 3'* / 7n Im 7w. 




which can be simplified to 


Oli = 


li 


li 


7m + Z^fc=l Ik Ik 


(40) 


(41) 


Thus, each resource partition fraction is equal to the fraction of the particular SNR over the total SNR. 


Substituting the above into (38) will always give the maximum sum rate 


N 


n = iog2 1 + X] 7 i 


(42) 


2=1 
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